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EVIDENCE FOR STELLAR SUBSTRUCTURE IN THE HALO AND OUTER DISK OF M31 
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ABSTRACT 

We report the discovery of significant stellar substructure in the halo and outer disk of our nearest 
large galactic neighbour, M31. Our deep panoramic survey with the Isaac Newton Telescope Wide Field 
Camera currently maps out an area of « 25 LJ° around M31, extending along the semi-major axis to 
55 kpc, and is the first to allow an uninterrupted study of the density and color distribution of individual 
red giant branch stars across a large fraction of the halo of an external spiral galaxy. We find evidence for 
both spatial density and metallicity (as inferred from colour information) variations, which are often, but 
not always, correlated. In addition to the previously reported giant stellar stream (Ibata et al. 2001b), 
the data reveal the presence of significant stellar overdensities at large radii close to the south-western 
major axis, in the proximity of the very luminous globular cluster Gl, and near the north-eastern major 
axis, coinciding with and extending beyond the previously-known 'northern spur'. The most prominent 
metallicity variations are found in the southern half of the halo, where two large structures with above 
average metallicites are apparent; one of these coincides with the giant stellar stream while the other 
corresponds to a much lower-level stellar enhancement. Our findings contrast with, but do not conflict 
with, past studies of the M31 halo and outer disk which have suggested a rather homogeneous stellar 
population at large radius: the bulk of our newly-detected substructure lies in the previously-uncharted 
far outer regions of the galaxy. We discuss the possible origin of the substructure observed and the 
implications it has for constraining the galaxy assembly process. 



Subject headings: galaxies: individual (M31)~ galaxies: evolution 
galaxies: stellar content - galaxies: structure 



galaxies: halos - Local Group 



1. INTRODUCTION 

Within the context of the hierarchical clustering theory 
for structure formation, large disk galaxies like the Milky 
Way and M31 arise from the merger and accretion of many 
smaller subsystems and from the smooth accretion of in- 
tergalactic gas (eg. White & Rees (1978); Steinmetz & 
Navarro (2002)). While the former process governs the 
growth of the dark matter halo (by far the dominant mass 
component in a disk galaxy) and much (if not all) of the 
stellar spheroid, angular momentum-conserving collapse of 
the latter is the currently-favoured mechanism for forming 
the thin disk. It has been argued that the relative im- 
portance of these processes - i.e. smooth versus discrete 
accretion - in the assembly of a galaxy is the primary 
factor which determines the final morphology (eg. Baugh, 
Cole, & Frenk (1996); Steinmetz & Navarro (2002)). 

Cosmological simulations of galaxy formation have be- 
come increasingly sophisticated in recent years and have 
now attained sufficient resolution to begin to address the 
internal structures of galaxies. A key result to emerge from 



state-of-the-art simulations within the popular cold dark 
matter (CDM) framework is that the dark matter subhalos 
which merge to form massive galaxies are much more re- 
silient than previously thought (Klypin et al. 1999; Moore 
et al. 1999). Although dynamical friction and galactic tidal 
forces continually act to disrupt subhalos once they fall 
within the potential of the massive system, the dense cen- 
tral cores of the satellites appear to survive as distinct 
entities for at least several orbital timescales (eg. Hayashi 
et al. (2002)). The total mass in undisrupted cores is ex- 
pected to be quite small (<J 10% of the total mass of the 
system) with the bulk of the mass distributed much more 
smoothly, however several hundred of them are expected 
to reside within the virial radius of a typical galaxy like 
the Milky Way. 

Signatures of the hierarchical nature of galaxy assem- 
bly are expected to be most obvious in the properties of 
the halo stellar populations. At least some of the accreted 
satellites will have experienced significant star formation 
prior to and during incorporation into the final system. 
As these systems will span a range of masses and have ex- 
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perienced different dynamical histories within the massive 
galaxy potential, it seems highly likely that they will also 
be characterised by distinct star formation and chemical 
enrichment histories, the details of which will be imprinted 
on the ages and metallicities of their constituent stars. In- 
deed, studies of Milky Way satellite stellar populations 
reveal a surprising variety of complex star formation his- 
tories, with no two dwarfs sharing a similar evolution (Ma- 
teo 1998). During the accretion of a satellite galaxy onto a 
massive host, tidal forces first shock the dark matter, then 
the stars, of the satellite making them gravitationally un- 
bound. Tidal debris will be deposited in both a leading 
and a trailing stream which, depending on the satellite's 
orbit, initial mass and sphericity of the host potential, can 
maintain spatial and kinematic coherence within the halo 
for many Gyr (Johnston, Hernquist, & Boltc 1996; Helmi 
& White 1999; White & Springel 2000; Bullock, Kravtsov, 
& Weinberg 2001). It is therefore expected that the stellar 
halos of galaxies which form hierarchically should posses 
significant spatial and metallicity substructure in the form 
of disrupted satellites (eg. dwarf galaxies, globular clus- 
ters) and their stellar detritus. This structure will be par- 
ticularly apparent in the outer regions of the halo, where 
the dynamical mixing timescales are long. Furthermore, 
as there is no reason to believe that the accretion history 
will be exactly the same in different galaxies (for example, 
the number and masses of accreted satellites and their mix 
of stellar populations), one might expect to sec noticeable 
differences in the properties of galactic stellar halos, even 
between systems of similar morphological type. 

It has also been postulated that some fraction of the ac- 
creted satellites lack a luminous stellar component, either 
due to feedback from early star formation which expels the 
gas before significant star formation can occur (eg. White 
& Rees (1978); Dekel & Silk (1986)) or from the inhibi- 
tion of gas cooling due to an ionizing UV background (eg. 
Bullock, Kravtsov, & Weinberg (2000)). The main rea- 
son for invoking these scenarios is the striking discrepancy 
between the numbers of satellites predicted to lie within 
the dark halo of a massive disk galaxy and the numbers 
of dwarf companions actually observed around the Milky 
Way and M31 (Klypin et al. 1999; Moore et al. 1999). If 
correct, detecting these CDM substructures through stel- 
lar density enhancements becomes very difficult, if not im- 
possible. On the other hand, such clumpy structures will 
still exert an influence on the the host galaxy by caus- 
ing tidal heating and/or distortion of the thin stellar disk, 
particularly in the fragile outer regions (Moore et al. 1999; 
Font et al. 2001). 

Searching for and studying stellar substructure in the 
outskirts of disk galaxies provides an important test of 
CDM models of structure formation. Studies of the Milky 
Way stellar halo have so far provided encouraging sup- 
port for hierarchical galaxy assembly Evidence includes 
the discovery of the tidally-disrupted Sagittarius dwarf 
galaxy and associated stellar streams (eg. Ibata, Gilmore, 
& Irwin (1995); Mateo, Olszewski, & Morrison (1998); 
Majewski et al. (1999); Yanny et al. (2000); Ibata et al. 
(2001a,c)), the phase-space clumping of halo stars in the 
solar neighbourhood (Helmi et al. 1999) and the tidal tails 
emanating from halo globular clusters (eg. Grillmair et 
al. (1995); Odenkirchen et al. (2001); Leon, Meylan, & 



Combes (2000)) and dwarf Spheroidal galaxies (eg. Irwin 
& Hatzidimitriou (1995); Majewski et al. (2000)). Unfor- 
tunately, our location within the disk sometimes renders 
interpretation of structures revealed through star counts 
rather difficult (eg. Newberg et al. 2002) and underscores 
the need for an extragalactic perspective. It is thus clearly 
desirable to test whether the halos of other disk galaxies 
also exhibit substructure and to investigate how the nature 
of this substructure (as well the nature of the field popula- 
tion resulting from totally disrupted satellites) varies with 
the properties of the host galaxy. 

Searches for low surface brightness tidal and/or extra- 
planar features have been carried out around several ex- 
ternal edge-on disk systems but to date only one unam- 
biguous detection has been made, in NGC 5907 (Shang et 
al. 1998). On the other hand, such studies are technically 
very difficult, and inaccurate flat-fielding and/or low-level 
scattered light can easily mask real signal at these levels 
(Sy ;> 28 magnitudes arcsec~ 2 ). A more preferable search 
technique involves using resolved star counts to map out 
the spatial density structure in nearby galaxy halos, how- 
ever, with ground-based telescopes, this method is cur- 
rently limited to galaxies within and around the Local 
Group. As these galaxies also subtend the largest angu- 
lar sizes on the sky, quantitative study requires wide-field 
CCD mosaic cameras as well as sophisticated processing 
techiques capable of dealing in an optimal way with large 
amounts of data, both of which have only become available 
in recent years. 

Using the Wide-Field Camera on the Issac Newton Tele- 
scope (INT WFC), we are carrying out a panoramic imag- 
ing survey of our nearest large neighbour, M31. The con- 
tiguous nature of our survey allows us to distinguish local 
density enhancements in M31's halo and outer disk from 
fluctuations in background galaxy counts and the fore- 
ground distribution of Galactic stars. The present paper 
reports interim results from this survey, namely the discov- 
ery of significant spatial density and metallicity substruc- 
ture in the red giant population of the halo and outer disk. 
The stellar halo of M31 has long been known to have rather 
different properties from the Milky Way halo, despite the 
overall similarity in the global properties of the galaxies. 
In particular, ground-based and HST studies indicate the 
M31 halo is roughly an order of magnitude more metal- 
rich (eg. Morris et al. (1994); Rich et al. (1996); Holland, 
Fahlman, & Richer (1996); Durrell, Harris, & Pritchet 
(2001); Reitzel & Guhathakurta (2002)) and denser and/or 
larger (Reitzel, Guhathakurta, & Gould 1998) than that of 
the Milky Way. Our INT WFC survey reveals significant 
stellar substructure in the outskirts of M31 and provides 
the tantalizing suggestion that these halo differences may 
be due to a more active accretion/merger history within 
the M31 subsystem. Preliminary results, drawn from the 
first phase of our survey, were reported in Ibata et al. 
(2001b). 

2. OBSERVATIONS 

The Wide Field Camera on the 2.5m Isaac Newton Tele- 
scope is a 4-chip EEV 4kx2k CCD mosaic camera which 
images «0.29 square degrees per exposure (Walton et al. 
2001). On the nights of 3-9 of September 2000, 9-16 Oc- 
tober 2001 and 13 November 2001, we used this camera 
to image 91 contiguous fields (corresponding to « 25 D°) 
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in the outer disk and halo of M31. Coverage currently ex- 
tends to 4.0° (« 54 kpc) and 2.5° (« 34 kpc in projection) 
along the major and minor axes respectively. 

Images were taken in the equivalent of Johnson V and 
Gunn i bands under mainly good atmospheric conditions, 
with 85% of the fields taken in photometric conditions 
with seeing better than 1.2". The exposure time of 800- 
1000s per passband per field allowed us to reach i = 23.5, 
V= 24.5 (S/N«5) and is sufficient to detect individual red 
giant branch (RGB) stars to My ~ and main sequence 
stars to My w — 1 at the distance of M31. Several fields 
taken in poorer conditions were reobserved and coadded as 
necessary to give an approximately uniform overall survey 
depth. 

All of the on-target data plus calibration frames were 
processed using the standard INT Wide Field Survey 
pipeline provided by the Cambridge Astronomical Sur- 
vey Unit (Irwin & Lewis 2001). This package provides 
the usual facilities for instrumental signature removal, in- 
cluding in this case dcfringing of the i-band data, plus 
tools for object catalogue creation, astrometric and photo- 
metric calibration, morphological classification and cross- 
matching catalogues from different observations. The 
pipeline processing provides internal cross-calibration for 
the four CCDs at a level better than 1% within each point- 
ing. Field-to-field variations in photometric zero-points 
were calibrated and cross-checked using a combination 
of multiple nightly photometric standard sequence obser- 
vations and the overlap regions between adjacent WFC 
pointings. The overall derived photometric zero-points for 
the whole survey are good to the level of ±2% in both 
bands. The full details of the survey strategy, data pro- 
cessing and calibration will be presented elsewhere (Irwin 
et al. , in prep). 

Objects were classified as noise artifacts, galaxies, or 
stars according to their morphological structure on all the 
images. As a sanity check, we compare our extended 
source counts in the outer halo fields with the deep V 
and I band galaxy counts of Smail et al. (1995). Over the 
magnitude ranges 22.5 <V< 23.5 and 20.5 < i < 21.5, 
we detect m 11500 and w 7300 galaxies per square degree 
respectively, which can be compared with the « 10000 per 
square degree in each band predicted using Smail et al. 
(1995) and indicates that we are not mis-classifying large 
numbers of faint galaxies. In the outer halo fields, we typi- 
cally detect equal numbers of stars and galaxies within the 
magnitude and color ranges of interest (see Section 3.1 
and Figure 1). Of the extended sources, approximately 
20% are 'compact' in the sense of being within the 3-5 
sigma range of the stellar boundary in the classification 
statistic and having an ellipticity of < 0.4. Making the 
plausible assumption that half of these are genuine galax- 
ies (the other half being genuinely stellar) leads us to we 
expect that the contamination due to mis-classified barely 
resolved field galaxies is small and generally considerably 
less than 10% of the total number of detected sources. In- 
deed, the overdense regions we will discuss in this paper 
have stellar densities at least a factor of two higher, and 
the average star counts are a factor of three higher, than 
that in the outer halo fields. 

Stellar contamination from the Galactic foreground 
steadily increases to the northeast, due to the proxim- 



ity of the Galactic Plane, rising smoothly from an aver- 
age contamination of « 13000 stars per square degree at 
the southwest extremity of the survey to «20000 stars per 
square degree at the northeast extremity (integrated over 
all magnitudes). This foreground variation, coupled with 
our current lack of suitable comparison fields uncontami- 
nated by M31, constrains to some extent detailed quanti- 
tative analysis of the spatial and metallicity distribution 
of the outer parts of the halo. 

3. RESULTS 

3.1. Spatial Density Variations 

Our INT WFC survey of M31 provides the first oppor- 
tunity to make an uninterrupted study of the properties 
of the resolved stars across a large fraction of an external 
disk galaxy. 'Stellar' sources in our catalogue consist of 
M31 red giants, M31 upper main sequence stars, Galactic 
foreground stars and unresolved background galaxies. In 
order to examine the spatial distribution of various stel- 
lar populations within M31, we apply a series of magni- 
tude and color cuts which are designed to isolate stars in 
different regions of the color-magnitude diagram (CMD). 
Figure 1 illustrates the main selection criteria adopted in 
our analysis. We define the 'blue RGB' as stars with 
20.5 < i < 22.5; 8.7 - 0.4 i < V - i < 23.5 - i and the 'red 
RGB' as stars with 20.75 < i < 21.75; V - i > 2.0. These 
magnitude and color ranges are consistent with those ex- 
pected for red giant stars at the distance of M31. We also 
overlay several Galactic globular cluster fiducials spanning 
a range in metallicity. Comparison between our selection 
criteria and these fiducials demonstates that our RGB cuts 
isolate the metal-poor (—2.0 <^[Fe/H]<j —0.6) and metal- 
rich (—0.6 <^[Fe/H]<; —0.3) sections of the giant branch 
respectively. 

Figure 2 shows the standard coordinate projection of 
the surface density distribution of blue and red RGB stars 
across our current w 25 D° survey area. Each detected 
source classified as stellar on the z-band and falling within 
the aforementioned magnitude and color limits is encoded 
as a point on these diagrams. The lower magnitude limits 
are conservatively set at roughly one magnitude brighter 
than the survey 5-cr detection threshold to mitigate the ef- 
fects of varying completeness. The immediate impression 
from Figure 2 is the striking non-uniformity of the stellar 
distribution at large radii. The large-scale morphology of 
stars in the blue RGB map is that of a significantly flat- 
tened inner halo structure, while the bulk of the stars in 
the red RGB map appear associated with the outskirts 
of the stellar disk. Numerous stellar enhancements are 
present in both maps, with some features appearing more 
conspicuous in one color cut than another. The giant stel- 
lar stream discovered during the first phase of our survey 
(Ibata et al. 2001b) is apparent as an enhancement close 
to, but distinct from, the southern minor axis and extend- 
ing out to ~ 40 kpc. While visible in both maps, the 
feature appears considerably sharper in the distribution of 
red RGB stars. In addition, significant stellar overdensi- 
ties are seen at large radii lying close to both major axes 
(at 1.5° W, 1.8° S and 0.8° E, 1.8° N). A lower intensity 
enhancement is visible in the form of a diffuse extension 
northeast of the center (towards 1.7° E, 0.7° N). Very faint 
structure, some of which is at the limit of detectability, is 
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seen emanating from the northern side of the disk. No 
correction has been applied to the maps for Galactic fore- 
ground contamination, however we have verified that the 
distribution of stars with properties expected to belong to 
this population (selected from the box indicated in Figure 
1) is smooth across the entire survey area. It is there- 
fore unlikely that Galactic foreground variations are caus- 
ing any of the structure seen in Figure 2. Likewise, the 
Galactic extinction towards the outskirts of M31 is rela- 
tively uniform, ranging from E(V— i)=0.06-0.11 (Schlegel, 
Finkbciner, & Davis 1998). There is no evidence for sig- 
nificant clumping in the distribution of extended sources 
(i.e. background galaxies) across our survey area and, in 
any case, galaxy clusters and other large scale structure 
would be unlikely to dominate the counts at the magni- 
tudes, colors and angular scales of the observed features. 
Indeed, Couch ct al. (1993) show that the galaxy-galaxy 
correlation function is very low beyond 0.1° at the magni- 
tudes of relevance for our survey. We thus conclude that 
the spatial density enhancements reflect real substructure 
associated with M31's halo and outer disk. 

The substructure located near the southwestern major 
axis is of particular interest because it is located in the 
proximity of the very luminous globular cluster Gl (Mey- 
lan et al. 2001). This clump, which we will hereafter refer 
to as the 'Gl clump' although the nature of the association 
with Gl, if any, is presently unclear, is located at ~ 35 kpc 
in projected radius and has a physical size of <~ 0.7° x 0.5°, 
or 10 x 7 kpc at the distance of M31. Compared with the 
mean density of RGB stars at similar radii around M31, 
we calculate that the Gl clump is overdense by a factor of 
«4. Integrating the excess RGB population over this re- 
gion relative to nearby comparison fields and assuming a 
luminosity function similar to the field population studied 
by (Rich et al. 1996) gives estimates for the total apparent 
magnitude of this feature of my — 12.1 and mi — 10.5, 
with average surface brightnesses of around 28.5 and 27.0 
magnitudes arcsec -2 respectively. At an average extinc- 
tion of E(B-V) = 0.07 this is equivalent to absolute mag- 
nitudes of My = -12.6 and M, t = -14.1. While the mag- 
nitude and color of the clump are comparable, within the 
errors, to those of dwarf satellites in the M31 and Milky 
Way sub-groups, the measured surface brightness is sev- 
eral magnitudes fainter than any of the currently known 
systems; on the other hand, it is similar to that measured 
for the low surface brightness ring in NGC 5907 (Shang et 
al. 1998). 

The density enhancement lying near the northeastern 
axis coincides with and extends beyond the 'northern 
spur', first remarked upon by Walterbos & Kennicutt 
(1988) (hereafter WK88) as faint light bending away from 
the major axis (see also the deep images of Innanen et al. 
(1982)). WK88 questioned whether this feature was actu- 
ally associated with M31, or due to a Galactic reflection 
nebula: our color-magnitude diagrams of this region un- 
ambiguously place the excess stars at the same distance 
as M31. The feature extends to a projected radius of 
~ 30 kpc at angles up to ~ 20° off the major axis. The 
northern spur region is about a factor of 1.5 - 2 times 
more overdense than the Gl clump and thus has a mean 
surface brightness w 0.5 — 0.8 magnitudes arcsec -2 higher. 
The northern spur region also appears as a low level en- 



hancement on a spatial density plot of stars selected to lie 
just above and redward of the RGB tip (see Figures 1 and 
3). Such stars are likely to be intermediate-age, moderate- 
metallicity ([Fe/H] > -0.7; ages ^3-8 Gyr) thermally- 
pulsing asymptotic giant branch (AGB) stars (Girardi et 
al. 2000), and their presence is taken as evidence for an 
extended epoch of star formation in this region. There is 
also a very tentative detection of AGB stars in the vicinity 
of the Gl clump, however all other substructure appears 
devoid of this younger population. While the effects of 
crowding will lead to some spurious structure on the AGB 
map (eg. the strong detection of NGC 205), this is unlikely 
to be important in the diffuse outer regions of the galaxy 
under study here. 

Perhaps surprisingly, our current map of the RGB den- 
sity around M31 reveals no obvious northern counterpart 
to the giant stellar stream discovered near the southern 
minor axis (Ibata et al. 2001b). We will return to this 
issue in the discussion below. 

3.2. Metallicity Variations 

Further insight on the nature of the stellar density vari- 
ations in M31's halo and outer disk comes from analysis 
of the spatial variation of mean RGB color. Figure 2 al- 
ready reveals that the morphology of the M31 substructure 
has a color dependence. This is further illustrated in Fig- 
ure 4, which shows color-magnitude diagrams for fields in 
the giant stellar stream (right) and the Gl clump (left). 
Comparison with Galactic globular cluster fiducials reveals 
that the mean RGB color varies between these pointings, 
being slightly more metal-rich than 47 Tuc in the stel- 
lar stream and approximately equal to 47 Tuc in the Gl 
clump. The width of the Gl clump RGB also appears 
somewhat narrower than that of the stellar stream, indi- 
cating a small intrinsic metallicity and/or age dispersion. 
In order to quantify the RGB color variation as a function 
of position over the survey area, it is first necessary to 
make a statistical correction for the foreground Galactic 
contamination. Three outer halo fields were used to de- 
fine a reference foreground population (albeit with a small 
halo contamination still present in these fields) and the 
scale factor required to match the target field z-band lu- 
minosity function brighter than i = 20 computed using 
the foreground box in Figure 1. For each field, we then 
calculate the foreground-corrected median color distribu- 
tion in the magnitude range 21 < i < 22 projected in a 
coordinate system orthogonal to the locus of the ensemble 
RGB for M31. This color measure quantifies the shift in 
the mean RGB locus as a function of field pointing. Figure 
5 shows a color-coded plot of this variation. The 'average' 
halo RGB color is denoted here as green; assuming an 
old stellar population, this corresponds to a metallicity 
of [Fe/H]~ —0.7, similar to 47 Tuc. Progressively bluer 
shades indicate bluer mean RGB colors while progressively 
redder shades indicate redder mean RGB colors. 

We adopt the view here that metallicity variations are 
the dominant driver of these color variations. While age 
variations may contribute somewhat to the observed be- 
haviour, there is yet to emerge any compelling evidence 
for a luminous AGB component in M31's extended halo 
which would accompany a young-to-intermediate age pop- 
ulation (eg. Rich et al. (1996); Holland, Fahlman, & Richer 
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(1996); also see Figure 3). For reference, an age variation 
of ~ 15 Gyr at [Fe/H] = —0.7 would be required to ex- 
plain the entire color variation seen over the survey area 
(Girardi et al. 2000). Based on the assumption of a uni- 
formly old age for the M31 halo, comparison with Galactic 
globular cluster tracks and theoretical isochrones indicates 
that color index variations of ±0.2 in V— i (the full range 
of colors used in Figure 5) correspond to a metallicity vari- 
ations ranging from -0.5 dex below the mean to 0.25 dex 
above the mean. As the relationship between RGB color 
and metallicity is highly non-linear across the range of in- 
terest (see Figure 4), the quoted spread in metallicities is 
unfortunately rather uncertain. The observed color varia- 
tions are well in excess of our photometric errors at these 
magnitudes, and also significantly larger than the zero- 
point uncertainties. Although our strategy of bootstrap- 
ping the calibration of non-photometric pointings from ad- 
jacent fields could introduce some correlated color varia- 
tions, the size of this effect is also estimated to be rather 
small, amounting to no more than a few hundredths of a 
magnitude. 

Inspection of Figure 5 reveals that while there are no ob- 
vious radial metallicity gradients present in the halo, there 
are significant chemical inhomogeneities on large (i.e. sev- 
eral kpc) scales. Further, comparison of Figure 5 and Fig- 
ure 2 indicates that while substructure in the metallicity 
map often corresponds to overdensities on the spatial den- 
sity map, this is not universally true. The most prominent 
chemical substructure is found on the southern side of the 
galaxy, in the region of the stellar stream and further to the 
north on the same side, where only a diffuse extension of 
the halo is present on the spatial density map. Both these 
regions appear significantly more metal rich than the aver- 
age halo and it is tempting to speculate that the two might 
be related, possibly tracing out the projected orbit of the 
object producing the stellar stream. We note that line-of- 
sight depth through the halo may cause intrinsic metallic- 
ity variations to be somewhat 'washed out', depending on 
how localised the substructure is in space and how large 
the overdensity is. The true metallicity variations in the 
halo arc thus likely to be somewhat underestimated in Fig- 
ure 5. 

Another apparently metal-rich region is the northern 
spur; at first glance this might appear to be even more 
metal rich than the outer stellar disk (the edge of which 
is approximately delineated by the inner ellipse) however 
the disk fields contain OB stars which make them ap- 
pear bluer, leading to spuriously low metallicity estimates. 
Similarily, age effects also complicate the interpretation of 
the spur, where we have evidence for the presence of an 
intermediate-age population (Figure 3). The metal-rich 
nature of this feature appears rather robust though, given 
that the young population will cause the instrinsic stellar 
colors to appear bluer, not redder. 

The Gl clump does not appear particularly distinct on 
the metallicity map, contrary to its appearance on the 
spatial density map, with the RGB stars in this region ap- 
pearing just slightly more metal poor than the background 
halo. 



4. DISCUSSION 



4.1. Comparison to Previous Studies of the M31 Outer 
Disk and Halo 

4.1.1. Optical Studies 

Many previous studies have addressed the stellar popu- 
lation characteristics of the halo and outer disk of M31, 
however these have almost exclusively sampled only a 
few discrete locations or have taken a panoramic, but 
much shallower view. Early quantitative wide area survey 
work of M31 was based on photographic plates (eg. Hodge 
(1973); Innanen et al. (1982); Walterbos & Kennicutt 
(1988)). Although these studies clearly revealed, among 
other things, the anti-symmetric warping of the outer disk 
and the tidally-induced twisting of the outer isophotes of 
the dwarf companions (see also Choi, Guhathakurta, & 
Johnston (2001)), they did not go deep enough to directly 
resolve the stellar populations. 

More recent work has focused on detailed studies of re- 
solved stars in several (generally small) fields sampling 
the outer disk and halo (see for example Mould & Kris- 
tian (1986); Morris et al. (1994); Rich et al. (1996); Hol- 
land, Fahlman, & Richer (1996); Reitzel, Guhathakurta, 
& Gould (1998); Durrell, Harris, & Pritchet (2001); Fer- 
guson & Johnson (2001); Sarajedini & Van Duyne (2001); 
Reitzel & Guhathakurta (2002)). Both ground-based and 
HST studies have found a dominant population which 
is significantly more metal-rich than that of the Milky 
Way halo and which is characterised by a considerable 
intrinsic dispersion. Specifically, measurements indicate 
[Fe/H] ~ —0.7, comparable to 47 Tuc, with spread of 
nearly 2 dex for M31, compared to [Fe/H] <~ —1.5 for 
the Milky Way. Detailed metallicity distributions have 
now been calculated for stars in several fields, all of which 
are found to exhibit a similar shape, consisting of a peak 
at the metal-rich end and an extended tail towards lower 
metallicities (Holland, Fahlman, & Richer 1996; Durrell, 
Harris, & Pritchet 2001; Sarajedini & Van Duyne 2001; 
Reitzel & Guhathakurta 2002). The fraction of stars in 
the metal-poor component is estimated to be 25-50% with 
no compelling evidence for a strong metallicity gradient 
and/or inhomogeneities in the halo. Additionally, Reitzel, 
Guhathakurta, & Gould (1998) find that M31's halo is 
roughly an order of magnitude larger and/or denser than 
the Galactic halo, though this conclusion is drawn from 
star counts in only three fields. 

Analysis of individual color-magnitude diagrams from 
our WFC survey indicates a mean metallicity in broad 
agreement with the findings of previous studies (see Fig- 
ure 4; also Ibata et al. (2001b)) and reveals no obvious 
metallicity gradient in the halo out to at least 50 kpc. 
We do, however, find evidence for significant spatial and 
metallicity substructure in the halo. Figure 6 shows the 
location of published HST fields and a few recent ground- 
based fields overlayed on our survey area. It can be seen 
that these fields generally sample either the inner halo of 
M31 or at large radius along the southern minor axis, nei- 
ther of which exhibit significant substructure in our survey. 
The chemical and spatial homogeneity of the halo that has 
been inferred from these studies is therefore not surpris- 
ing. Most of the substructure our survey has uncovered 
lies in the previously- uncharted outer regions of M31. 

An exception to this is the study by Morris et al. (1994) 
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which resolved the stars in a small section of the northern 
spur and provided the first evidence that the stellar over- 
density in these parts consisted of the fairly characteristic 
metal-rich population seen elsewhere in the outskirts of 
M31. Like us, these authors also found evidence for a sig- 
nificant AGB component lying above the RGB tip. In ad- 
dition, Rich ct al. (1996) obtained an HST/WFPC2 CMD 
of the field population around the globular cluster Gl, cor- 
responding to a location near the very edge of the stellar 
overdensity discovered in our WFC survey, and found the 
RGB to be well-matched by the 47 Tuc fiducial. They also 
noted that the luminosity function of the field population 
could not be fitted by template globular cluster luminosity 
functions. While the results from these studies are in good 
qualitative agreement with those presented here, their lack 
of depth and/or areal coverage adds little further insight 
into the nature of the stellar populations associated with 
the substructure. 

4.1.2. HI Studies 

The most extensive published wide-field survey of HI in 
M31 is still that of Newton & Emerson (1977) who mapped 
regions along the major axes out to ~ 35 kpc. Apart from 
the warp in the disk at large radius, the outermost HI con- 
tour looks remarkably smooth and uniform with no sign of 
perturbation in the direction of any of the halo substruc- 
ture (sec their Figure 12a). Furthermore, while the optical 
and HI disks are warped in the same direction, the optical 
warp appears to begin at smaller radii than that of the HI 
and exhibit a greater deviation from the plane (see also In- 
nanen et al. (1982); Walterbos & Kennicutt (1988)). This 
holds true along both major axes, but seems particularly 
apparent in the north-east. 

An interesting high velocity HI cloud lying near the 
northern minor axis of M31 was discovered by Davies 
(1975). The feature lies approximately 1.5° (« 20 kpc) 
north-west of the center of M31 at (0.9°W, 1.1°N) , has 
an extent of ~ 1° and a heliocentric velocity of —450 km/s 
(see Figure 8). It lies just above NGC 205 and, rather cu- 
riously, along the projected extension of the giant stellar 
stream discovered in the southern half of M31. Based on 
the velocity of the cloud, it would appear unlikely that it 
has any connection to the north-west disk, nor for that 
matter to either M32 (Y heUo = -197 km/s) or NGC 205 
(Vft e /i = —242 km/s). The cloud has a peak column den- 
sity of 4 x 10 19 atom cm -2 , a gaseous velocity dispersion 
of 9 km/s and, if at the distance of M31, a gas mass of 
5 x 10 6 M Q ; it thus represents a fairly typical compact 
high velocity cloud. Inspection of CMDs and detailed spa- 
tial density and metallicity maps centered on the cloud 
position reveal nothing unusual at this location. An al- 
ternative explanation for the cloud is that it is very local 
and related to the Magellanic Stream (Davies 1975); the 
non-detection of an associated stellar counterpart to the 
cloud may be more consistent with this interpretation. 

Newton & Emerson (1977) also discovered an apparently 
detached HI cloud at large radius along the north-eastern 
major axis of M31. Located at a projected distance of 
~ 3° or ~ 34 kpc, this cloud has a velocity which differs 
by 100 km/s from that expected for local disk rotation. 
The global properties of the cloud are similar to those in- 
ferred for the Davies cloud. The cloud lies at (1.5°E, 2°N) 



from the center of M31 and near our current survey limit 
in this quadrant; our maps reveal no obvious spatial sub- 
structure here, but the mean metallicity of the RGB stars 
is somewhat above average (see Figure 5). The HI fea- 
ture may be a signature of a major disk disturbance in the 
north-east quadrant, possibly related to the northern spur 
feature which is located roughly 0.5° further in. 

4.2. The Origin of the Substructure 

Substructure in the form of accreted and/or disrupted 
satellites is a generic prediction of hierarchical galaxy as- 
sembly. Even if the satellites are accreted before significant 
star formation has occurred, they will still interact grav- 
itationally with the stellar disk of the host, causing tidal 
heating, distortion and possibly even warping. Our WFC 
survey has led to the discovery of significant spatial and 
chemical substructure in the outskirts of M31, the most 
prominent of which is characterised by Figure 7. We first 
address whether any of this substructure is directly asso- 
ciated with the luminous satellite companions of M31. 

M31 is known to possess approximately 15 satellites, of 
which the dwarf ellipticals (dE) M32 and NGC 205 are 
the most proximate, lying at projected distances of 5 and 
9 kpc respectively, and among the most luminous (Mateo 
1998). Figure 8 illustrates where the closest satellite com- 
panions lie in projection with respect to M31. Apart from 
the 'clump' of stars close to the luminous cluster Gl and 
the northern spur, the nature of which are not yet clear, 
our survey has failed to reveal any obvious (i.e. still in- 
tact) additional satellites out to ~ 3° from the center of 
M31. For reference, all known Milky Way satellites (in- 
cluding the Sagittarius dwarf) have surface brightnesses 
and angular extents that would have rendered them easily 
detectable in our survey, if lying within a projected radius 
of « 55 kpc of M31. 

M32 and NGC 205 are known to exhibit anomolous 
surface brightness profiles and/or outer isophote twists 
suggestive of tidal interaction and disruption (eg. Hodge 
(1973); Choi, Guhathakurta, & Johnston (2001), see also 
Figure 9). In addition, both display a variety of puzzling 
properties for galaxies of the dE class. M32 exhibits a 
strong intermediate-age stellar component, in addition to 
a classical old stellar component, and shows evidence for 
a large metallicity spread with a mean just below solar 
(Grillmair et al. 1996; Davidge et al. 2000). The signa- 
tures of an extended and complex star formation epoch 
are even more obvious in NGC 205, which possesses cold 
gas (both atomic and molecular) and dust as well as young 
stars (Lee 1996; Young & Lo 1997). Furthermore, the HI 
shows a well-defined velocity gradient (unlike the stars) 
and is confined to the very inner regions of the optical 
galaxy; this inconsistency between gas and stars suggests 
that the gas may have been recently captured, perhaps due 
to a passage through the disk of M31. The geometrical 
alignment of the two satellites with the southern stream, 
and the broadly similar metal abundances initially sug- 
gested that either or both of them could be responsible for 
the tidal feature (Ibata ct al. 2001b). With a significant 
fraction of the northern half of M31's halo now mapped, 
we have failed to identify any obvious continuation of the 
stream on the northern side of NGC 205. Given that tidal 
interactions give rise to both leading and trailing streams 
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from the disrupted satellite, such a feature would be ex- 
pected if NGC 205 were the origin of the stream. M32 
or a third system (either more distant, or now completely- 
cannibalized?) now appear the most promising candidates 
for the origin of the stream. We have previously remarked 
on the possible connection between the stream and the 
similarily metal-rich, but lower-level, stellar ovcrdcnsity 
near the north-eastern side of the disk (see Figure 7). If 
this association is correct, it implies the projected orbit 
of the disrupting satellite wraps tightly around the M31 
nucleus, and does not cross into the northern half of the 
halo of M31. On the other hand, the lack of a northern 
counterpart to the stream might simply reflect the fact 
that the current orbit of the satellite has yet to traverse 
into this part of the halo, or else that the stream in this 
region has already dispersed and is no longer visible as a 
spatial or chemical density enhancement. 

Other satellites of possible relevance for the interpre- 
tation of the M31 substructure are the dwarf Spheroidal 
systems, Andromeda I and Andromeda III. Both lie just 
beyond the boundary of our current survey to the south- 
west of M31, consistent with the general direction of the 
stream. The stellar populations in these systems are pre- 
dominantly metal-poor however, a fact that would seem 
to argue against them being the origin of the more metal- 
rich stream. Specifically, And I and And III have mean 
RGB mctallicities of [Fe/H] ~ —1.5 and —1.9 respectively 
(Da Costa et al. 1996; Da Costa, Armandroff, & Caldwell 
2002) compared to the value of [Fe/H] ^ -0.7 inferred 
for the stream. M31's largest satellite, M33, lies at a pro- 
jected distance of ~ 15° towards the south-east. There 
is no evidence at present to connect M33 with any of the 
substructure seen around M31, however the potential role 
this system could play within the M31 subsystem should 
be borne in mind. 

The extreme compactness of M32, coupled with the very 
high central surface brightness and the dwarf-like luminos- 
ity, has led many authors to suggest that M32 is the left- 
over core of a larger galaxy, evolved to its current state as a 
result of prolonged tidal 'harrassmcnt' within M31's halo 
(eg. Faber (1973); Bekki et al. (2001); Graham (2002)). 
Such a hypothesis may also explain the puzzling obser- 
vation that M32 has no globular clusters; a system this 
luminous would be expected to possess at least ~ 15 — 20 
of them (Harris 1991). The broad agreement between the 
stellar metallicity and dispersion in M32 and in the M31 
halo (Grillmair et al. 1996; Durrell, Harris, & Pritchet 
2001) leads one to speculate whether the former could be 
the origin of not only the stellar stream, but also much of 
the field halo. If the M31 halo was significantly polluted 
by stars from a fairly massive companion, this could also 
account for the higher density/size of the stellar halo (Re- 
itzel, Guhathakurta, & Gould 1998) compared with that 
of the Milky Way. For such a scenario to be viable, M32 
must have already made several revolutions around M31 in 
order to account for the high covering factor of the metal- 
rich population. Dynamical friction due to stars and dark 
matter in M31's halo will cause the orbit of M32 to spiral 
inwards with time, while any non-sphericity in the poten- 
tial and/or close passages to the disk will cause the orbit to 
precess. The significant velocity dispersion of M32 stars 
(van der Marel et al. 1994) implies that once pulled off, 



these stars will merge rather quickly into the halo popu- 
lation and only the most recently stripped stars will still 
appear as a coherent structure. Penarrubia, Kroupa, & 
Boily (2002) have shown that tidal debris from a disrupt- 
ing satellite will spread out even faster if the host halo is 
significantly flattened, with a timescale of 3-4 Gyr before 
total disruption of the companion. Detailed dynamical 
modelling of satellite orbits within the potential of M31 
will be presented in a future paper. 

Few constraints exist on the origin of the major-axis 
substructure at present. It would appear unlikely that 
the Gl clump and the northern spur are related to each 
other (for example, resulting from an object orbiting close 
to the plane of the disk) due to the different mean col- 
ors of the RGB stars in these regions and the presence of 
a prominent intermediate-age component in the spur but 
not in the clump. The Gl clump is particularly intrigu- 
ing given its proximity to the anomolous globular cluster 
Gl. Gl is among the most luminous and most massive 
globular clusters in M31 and is rather unique in having an 
intrinsic metallicity spread, suggestive of self-enrichment 
(Meylan et al. 2001). The outer isophotes of the cluster 
are distinctly elongated, and the direction of this elon- 
gation appears to be in the same sense as that of the 
Gl clump. Meylan et al. have argued that Gl might be 
the core of a disrupted dE and such a scenario is sup- 
ported by the fact that the metallicity reported for Gl 
([Fe/H]= —0.95) agrees quite well with that inferred for 
the stellar overdensity, i.e. slightly below that of the aver- 
age halo [Fe/H]~ —0.7. On the other hand, the significant 
spatial offset (0.5° or 7 kpc if at the distance of M31) 
seen between Gl and the peak of the stellar overdensity is 
difficult to understand in this picture, as is the fact that 
the luminosity and mean metallicity of the clump place 
it well off the luminosity-mctallicity relation defined by 
Local Group dE and dSph galaxies (eg. Mateo (1998)). 
Given the metallicity of the Gl clump as determined from 
mean RGB colour, the luminosity- metallicity relation pre- 
dicts a luminosity more than two magnitudes brighter than 
that observed. An alternative explanation for the feature 
is that it is a highly warped section of the far outer disk, 
perhaps perturbed or 'torn off' by a previous tidal interac- 
tion. The clump lies just beyond the outermost HI contour 
of Newton & Emerson (1977) which warps southward in 
this quadrant, compared to the northern offset displayed 
by the clump. If the disk interpretation is correct, the 
fact that the gas and stars currently exhibit different be- 
haviours here could imply the interaction happened a long 
time ago. Future photometric and kinematic observations 
will aid significantly in discriminating between these two 
possibilities. 

The location of the northern spur near the north-eastern 
major axis of M31 and in the same direction of the gaseous 
warp provides strong support for the association of this 
feature with a severe warp in the outer stellar disk. This 
interpretation would naturally explain the metal-rich na- 
ture of the population in the spur, as well as the presence 
of a young-to-intermediate age component, as these are 
properties which are also known to characterise stars in 
the outer parts of the disk (Morris et al. 1994; Ferguson & 
Johnson 2001)). Walterbos & Kennicutt (1988) noted that 
if the northern spur was an extension of the disk, it would 
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be the most extreme example of a warped stellar disk ever 
found. It would also be highly asymmetric. Tidal forces 
due to an infalling satellite could possibly excite such a 
strong warp (Binney 1992; Huang & Carlberg 1997), how- 
ever the identity and location of the putative perturber is 
not obvious at present. 

5. SUMMARY 

Our INT WFC survey currently maps an area of w 25 D° 
around M31, extending to a semi-major axis of 55 kpc, 
and allows the first uninterrupted study of the density 
and color distribution of individual red giant stars across 
a large fraction of an external disk galaxy. We have found 
evidence for both spatial and metallicity (as inferred from 
color information) substructure, which are often, but not 
always, correlated. In addition to the giant stellar stream 
reported by Ibata et al. (2001b), the data reveal the pres- 
ence of significant stellar overdensities at large radii close 
to the south-western major axis, in the proximity of the 
very luminous globular cluster Gl, and near the north- 
eastern major axis, coinciding with and extending beyond 
the 'northern spur'. The most prominent metallicity vari- 
ations are found in the southern half of the halo, where 
two structures with above average metallicities are appar- 
ent. One of these coincides with the aforementioned giant 
stellar stream, whereas the other corresponds with a much 
lower-level stellar overdensity. 

Our results contrast with, but do not conflict with, the 
findings of previous M31 stellar population studies that 
have found a rather homogeneous population in the halo 
and outer disk. These studies have mostly focused on inner 
halo fields or fields at large radius along the southern mi- 
nor axis, neither of which are found to exhibit significant 
substructure within the limits of our current survey. Of 
the newly-detected substructure, only the northern spur 
appears to have a counterpart in published HI maps (New- 
ton & Emerson 1977). Further, our data reveal no spatial 
or chemical overdensity coincident with the high velocity 
cloud discovered by Davies (1975) which lies just north 
of the satellite companion NGC 205, along the projected 
northern extension of the giant stellar stream. 

We have argued that the stellar substructure seen in 
the outskirts of M31 is likely due to a variety of past 
and ongoing interaction and accretion events. The lack of 
any obvious northern counterpart to the southern stellar 
stream suggests that either the disrupting satellite has yet 
to traverse into the northern half of the halo, or that the 
projected orbit wraps tightly around the center of M31. 
Indeed, there is tentative evidence to connect the stellar 
stream to the stellar and chemical enhancement seen to- 
wards the north end of the disk. Several pieces of evidence 
suggest that M32 could be the origin of the stream; these 
include the geometrical alignment of the two objects, the 
similarity in the mean metallicity and dispersion of their 
constituent stars and the wide range of peculiar properties 
exhibited by M32, all of which are consistent with it hav- 
ing once been a considerably more massive and luminous 
galaxy. One may even speculate that stars stripped from 
M32 have polluted a large fraction of the M31 halo, pro- 
viding a viable explanation for the long-standing puzzle of 
why M31's field halo is more metal-rich and denser/larger 
than that of the Milky Way. Alternatively, the stream 
may be due to another satellite, perhaps now completely 



cannabilized or else lying beyond our current survey limits. 

The stellar substructure located close to the major axes 
- the northern spur and the Gl clump - could plausibly 
result from tidal distortion and disruption of the outer 
disk due to, for example, the close passage of a massive 
satellite. If this explanation is correct, the interaction re- 
sponsible for producing the Gl clump must have occurred 
some time ago, given the lack of evidence for a strong 
intermediate-age component and the slightly lower than 
average metallicity in this structure. On the other hand, 
the remarkable proximity and overall alignment between 
the clump and very luminous globular cluster, Gl, sup- 
ports the alternative hypothesis that the two have a com- 
mon origin. Meylan et al. (2001) have proposed that Gl 
may be the core of a stripped dwarf elliptical galaxy; this 
hypothesis clearly warrants further consideration. How- 
ever, a satisfactory model must also account for the ~ 0.5° 
projected offset seen between the peak of the stellar over- 
density and the globular cluster. 

The Milky Way and M31 are the only two massive galax- 
ies for which we have detailed information about how the 
properties of individual stars vary across a large fraction 
of the halo. The mean metallicity and size/density of the 
stellar halos in these systems differ substantially and this 
has sometimes been cited as evidence for different forma- 
tion mechanisms (eg. (Durrell, Harris, & Pritchet 2001)). 
Our INT WFC survey reveals that, like the Milky Way, 
the stellar populations at large radii in M31 exhibit sig- 
nificant substructure. This substructure is expected in 
hierarchical CDM models of galaxy formation, cither as 
a direct result of disrupting satellites or as a by-product 
of tidal interactions between infalling satellites and the 
fragile outer disk. The fact that luminous, and clearly dis- 
torted, satellites lie close to the center of M31 suggests 
that the stellar halo differences between that galaxy and 
the Milky Way may be the result of a more active ac- 
cretion/merger history. Essential to the interpretation of 
our newly-detected substructure will be deep photomet- 
ric and spectroscopic observations of stars in the outskirts 
of M31. An approved Cycle 11 HST program with the 
newly-installed Advanced Camera for Surveys is aimed at 
obtaining deep color-magnitude diagrams, reaching sev- 
eral magnitudes below the horizontal branch, in several 
fields where we have detected halo substructure, and in 
two locations in the far outer disk. These data will provide 
more quantitative information about the ages and metal- 
licities of stars in the various substructures and will help 
to constrain the accretion/merger history of the galaxy. 
Spectroscopic observations with 8-m telescopes will be es- 
sential for providing kincmatical constraints to be used in 
modelling the satellite interactions, and for probing the 
shape and potential of M31's massive dark halo. 

Finally, it is worth noting that M31 is being used as 
a target in several large experiments that search for mi- 
crolensing events due to Massive Compact Halo Objects 
(MACHOs). Without further detailed modelling, it is im- 
possible to determine where the giant stellar stream lies 
with respect to center of M31. If it lies in front of the 
galaxy, it will produce a population of lenses that are more 
probable than the average M31 bulge star to lens M31 
bulge stars, whereas if it lies behind the galaxy, it will pro- 
duce an extra population of sources which are more likely 



to be lensed by M31 bulge stars. A careful appraisal of 
the stream, as well as other M31 stellar substructure, will 
therefore be required in order to interpret the microlensing 
event rates towards M31. 

We thank the Institute of Astronomy, Cambridge and 
the Anglo-Australian Observatory for hospitality during 
various collaborative visits. 
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Fig. 1. — An illustration of the selection criteria adopted to isolate stars in different regions of the color-magnitude diagram. The blue 
RGB (20.5 < i < 22.5; 8.7 - 0.4 i < V - i < 23.5 - i) and red RGB (20.75 < i < 21.75; V - i > 2.0) cuts are designed to select stars 
from the metal-poor and metal-rich sides of the red giant branch respectively. The AGB cut (19.75 < i < 21.00; V — i > 2.5) is designed to 
select intermcdiatc-age, moderate metallicity stars in the thermally-pulsing regime of the asymptotic giant branch. Ficld-to-ficld corrections 
for the varying Galactic foreground population (see Section 3.2) are calculated using star counts in the region defined by 18.0 < i < 19.5 
and 0.5 < V — i < 2.5. Also overlaid are the giant branch fiducial sequences of several Galactic globular clusters: from left to right these 
correspond to NGC 6397 ([Fe/H]= -1.9), NGC 1851 ([Fe/H]= -1.3), 47 Tuc ([Fc/H]= -0.7) and NGC 6553 ([Fe/H]= -0.3) taken from Da 
Costa & Armandroff (1990) and Sagar ct al. (1999) (see text). In placing the fiducials on the (V,i) plane, we have made use of the color 
equations derived for the WFS data (see http://www.ast.cam.ac.uk/~wfcsur/photom.html). 
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Fig. 2. — (a) A standard coordinate projection of the surface density of blue RGB stars across our current zss 25d° survey area. The inner 
and outer ellipses are drawn assuming a position angle of 38.1°, derived from analysis of the light distribution on a scanned Palomar Sky 
Survey plate. The outer ellipse denotes a flattened ellipsoid (aspect ratio 3:5) of semi-major axis length 55 kpc and indicates the current 
spatial extent of the survey. The inner ellipse has a semi-major axis of 2°(sd 27 kpc) and represents an inclined disk with i = 77.5; the optical 
disk of M31 lies well within this boundary. The few white blotches indicate regions contaminated by saturated stars. The dwarf companions 
M32 and NGC 205 lie at (0°,— 0.4°) and (—0.5°, 0.4°) respectively. Much substructure is seen at large radii, including the giant stellar stream 
and stellar overdensities at both extremes of the major axis. No corrections have been made for foreground or background contamination. 



13 



Fig. 2. — (b) Same as (a) except showing the surface density of red RGB stars. Note the lower Galactic foreground contamination on this 
map. Comparison with (a) clearly indicates that the morphology of the substructure varies as a function of color. 
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Fig. 3. — A standard coordinate projection of the surface density of stars lying redward of and above the RGB tip. These stars are likely to 
be intermediate-age, moderate metallicity thermally-pulsing asymptotic giant branch stars. Of the substructure visible in the RGB map, only 
the northern spur is clearly detected on the AGB map; there is also marginal enhancement in the vicinity of the Gl clump. No corrections 
have been made for foreground or background contamination. 
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Fig. 4. — Examples of two color-magnitude diagrams from our WFC survey. The left panel represents a field in the giant stellar stream 
and, the right panel, a field in the Gf clump. Over-plotted are the same globular cluster fiducial sequences as shown in Figure f. Most of 
the stars more luminous than the RGB tip are foreground contaminants. Stellar density varies between the CMDs as a result of the different 
overdensities of the features. We argue the different mean colors and widths of the RGB primarily reflect intrinsic metallicity variations. 
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Fig. 5. — A map of the variation in mean RGB color across M31 (see text for details). Each WFC pointing is represented as a color- 
coded polygon. The 'average' halo RGB color is denoted as green; assuming an old stellar population, this corresponds to a metallicity of 
[Fe/H]~ —0.7, i.e. similar to 47 Tuc. Yellow, orange and red regions indicate progressively redder mean colors, and hence higher mctallicities. 
Turquoise and light blue regions indicate progressively bluer mean colors, and hence lower metallicitics. The full metallicity range spanned 
by the plot (red through light blue) is £b 0.7 dex. Dark blue regions represent pointings where the mean stellar density was too low to make 
a definitive measurement of stellar color. 
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Fig. 6. — A diagram showing the locations of field population studies in the outskirts of M31. The squares show fields targeted by the 
ground-based studies of Morris et al. (1994); Reitzcl, Guhathakurta, & Gould (1998) and Durrell, Harris, & Pritchct (2001). The circles 
represent the HST/WFPC2 studies of Holland, Fahlman, & Richer (1996); Rich et al. (1996); Sarajedini & Van Duyne (2001); Ferguson, 
Gallagher, & Wyse (2000) and Ferguson & Johnson (2001). Only the Morris et al. 'spur' field and the Rich et al. Gl field lie close to halo 
substructure discovered in the present survey. 
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Fig. 7. — A cartoon illustrating the most prominent spatial and chemical substructure discovered in our WFC survey. A possible projected 
orbit of the giant stellar stream is indicated; the stream may connect to the enhancement seen in the spatial and chemical maps of the northern 
half of the disk. The contours delineate the approximate extent and orientation of the northern spur and Gl clump; the location of the Gl 
globular cluster is also indicated. 
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Fig. 8. — The projected distribution of the closest satellite companions around M31. Also shown is M31's nearest massive companion, M33 
and the position of the neutral hydrogen cloud detected by Davies (1975). As previously, the outer ellipse denotes the approximate limit of 
our current WFC survey. 
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Fig. 9. — Isopleth maps of M31 from an APM scan of a 75 min exposure Palomar Schmidt IllaJ plate taken by Sydney van den Bergh in 
1970. The lowest isophote contoured is at a level equivalent to B=27 mag arcsec -1 . The right-hand panels show the distorted outer isophotes 
of M32 and NGC 205 after substracting an elliptically-avcraged M31 profile. Cepa & Beckman (1988) quote tidal radii for M32 and NGC 205 
of 0.84 kpc (0.06°) and 1.93 kpc (0.14°) respectively, which corresponds to the relatively unperturbed inner parts of the profiles. The location 
of the Davies (1975) HI cloud is also indicated on the main figure. In the southern half of the galaxy, the M31 outer isophotes are distorted in 
the direction of the stream. Although the lowest contour level is at the limit of relibability of the photographic data, this particular feature 
is present in other digitised wide area photographic data (eg. Waltcrbos & Kennicutt (1988)). 
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